Exact symmetry breaking ground states for quantum spin chains 
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We introduce a family of spin-1/2 quantum chains, and show that their exact ground states break 
the rotational and translational symmetries of the original Hamiltonian. We also show how one 
can use projection to construct a spin-3/2 quantum chain with nearest neighbor interaction, whose 
exact ground states break the rotational symmetry of the Hamiltonian. Correlation functions of 
both models are determined in closed form. Although we confine ourselves to examples, the method 
can easily be adapted to encompass more general models. 
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The question of exact ground state of a given quantum 
many body system, i.e. a quantum spin chain, has always 
been fascinating and of great importance to physicists in 
all fields, from condensed matter physics and statistical 
mechanics to mathematical physics and quantum field 
theory [1-3]. Besides their inherent interest as a the- 
oretical laboratory, where the phenomenon of quantum 
phase transition [4] can be studied in detail and in ana- 
lytical form, quantum spin chains provide an important 
link with classical statistical mechanics on two dimen- 
sional lattice models, the latter is a subject with a rich 
and fascinating variety of cooperative and critical phe- 
nomena affording insights into many properties of real 
physical, chemical and biological systems. 

As prototypes of correlated many-body systems, the 
Hilbcrt space dimension of quantum spin chains grows 
exponentially with system size, rendering any numerical 
treatment of large systems extremely difficult. Therefore 
any exact solution of a new quantum spin chain, is not 
only valuable for its own specific characteristics, but may 
provide us with a plethora of mathematical techniques 
for constructing a whole new series of exactly solvable 
models. 

In this letter, we introduce a new model and tech- 
nique for constructing exact ground states of quantum 
spin chains which spontaneously break the symmetries of 
the parent Hamiltonian. The method is inspired by the 
original work of Affleck, Kennedy, Lieb and Tasaki who 
first introduced it for constructing exact invariant ground 
states for certain quantum chains [5]. Their method, de- 
signed in the context of valence bond solids, led in the 
course of time to the formalism of finitely correlated or 
matrix product states [6]. The model we introduce is 
a spin 1/2 quantum system on a ring, with interaction 
range equal to 4, having both translational and rotational 
symmetry, nevertheless we show that its exact ground 
states break both these symmetries. We also find the ex- 
act maximum energy state and some of the exact excited 
states, a progress which has not been reported for any 
of the models constructed in the AKLT or the matrix 
product formalism [7-10]. Finally we use this construc- 



tion to introduce another quantum spin chain, namely a 
spin 3/2 quantum chain with nearest neighbor interac- 
tion, with symmetry breaking ground states. We should 
stress that although we do this for a specific model with 
rotational and translational symmetry, our method can 
be adapted for construction of a broader class of similar 
models. 

Let us first fix our notation: we use |-|-), |— ) and \S) := 
-^([-1-,— ) — [ — ,+)) to denote respectively the spin up, 
spin down and a singlet state of two spins. Subscripts 
indicate the sites of the lattice, thus IS'12,+3) indicates 
a singlet on sites 1 and 2 and a spin up state on site 3. 
Consider the state 
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shown in figure (1). We require that |(/)o) be the ex- 
act ground state of a spin chain Hamiltonian H = 
Yfk=i^k,k-\-i,k-\-2,k-h3, where we use periodic boundary 
conditions. To this end, we demand that the local Hamil- 
tonian h he a, positive operator and annihilate |(/)o). A 
local Hamiltonian h finds four adjacent spins only in one 
of the mixed states shown in figure (1). Thus, we should 
have 



tr [h{\S){S\(g)\+){+\(g)I)]=Q, 

h[\+)\S)\+)] = 0, 

tr [h{I<»\+)(+\®\S){S\)]=0. 



(2) 



Since \S) is a singlet, the state |+)|S')|-f), can be written 
as a linear combination of states with total spin and 

1. Therefore a projector P2 which projects the product 
of four adjacent spins on the the sector with total spin 

2, automatically satisfies equations (2). However there 
is one other projector which satisfies equations (2), and 
its annihilation of these states is quite nontrivial. This is 
the projector Pq := where |x) (or its embedding 
on sites 1,2,3 and 4) is a singlet of four spins constructed 
as follows: 



IX) |S'i4, ^'23) + |5'i3, S24) 

To show this we use the simple facts 
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FIG. 1. The three equations in (2) correspond, respectively 
from left to right, to the annihilation of these local mixed 
states by h. For ease of comparison with the analysis in the 
text we have labeled the spins in all boxes as 1, 2, 3 and 4. 
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where 1,2,3 and 4 are any four different indices. From 
this identity one can easily derive other identities sim- 
ply by permuting indices in an appropriate way, namely 
(534 1514, 5*23) ~ -2-|'S'i2), or (534 1513, 5*24) — ^l'S'12). 
These two identities, immediately imply that (534|x) = 
0. Similarly one finds {Si2\x) = 0, hence the first and the 
last equations of (2) are satisfied. For the second relation 
of equation (2) we note that 



(XI'S'23) = (5*14, 5231523) + (5i3 , 524|523) lj^{Sl4\, (5) 

hence 



(xl + l,'S'23,+4> = 2^'^14| 



= 0. 



(6) 



This proves the assertion that Pq annihilates all three 
kinds of local states. Thus, we find that the local Hamil- 
tonian can be of the form h = 2JP2 + Pq, where J is a 
positive but arbitrary coupling constant, the factor of 2 is 
for convenience and one of the couplings, the coefficient 
of Po; has been re-scaled to unity. The final form of the 
Hamiltonian in terms of spin operators can be found as 
follows. First we note that i'^'' = 2 © l'^ © 0^ where the 
exponents indicate the multiplicities and Pq is the pro- 
jector to one of the above two zeros. Let us denote the 
sum of spin-j projectors in the above decomposition by 
Vj. Then, from 5 • 5 = J'O' + 1)"^^: where 5 is the 
spin operator of four consecutive sites, we have 



I ^Vo + Vi+V2, 
5 • 5 = 2Pi + 6^2, 
(5 • Sf = AVi + 36V2, 



(7) 



from which we obtain 
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^0 = 1- ^5-5 + -(5-5)2, 
^2 = -^5. 5 + 1(5. 5)2. 



(8) 



This readily gives P2 = {P2) in terms of spin operators. 
However, the individual projector Pq can not be obtained 
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FIG. 2. The states \<po), and \<p2)- 

in this way. Direct calculation from (3) is notoriously 
tedious. We can circumvent this problem by calculating 
the other projector Pq ~ |x^)(x^l and subtract it from 
Vo to obtain Pq. It is straightforward to check from (4) 
that the other singlet which should be perpendicular to 
Ix) is the following 
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This readily gives 



Po=Vo- \X^){X^) = T'o - (J - Si • S2)(i - S3 • S4). 

Putting all these together, and discarding overal con- 
stants, we find 
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-ff = ^ Si • Si+i(l - 2Sj+2 • Si+3) 



-iils,.S, + i±^(S,.S,)2, (10) 

where Si is the spin operator of one site and S; is the spin 
operator of a block of four spins at sites i,i + l,i + 2,i + 3. 
Note that the relative strength of consecutive couplings 
are approximately 1: 0.48: 0.06. 

The state \(j)o) is the exact ground state of H which 
breaks translational and rotational symmetry of the 
Hamiltonian. Clearly, it is not invariant under rotation, 
it is the top state of the total spin--y sector, other states 
which are degenerate with it in energy, are obtained by 
acting with total L_ — X)fc=i '^k operator, which does 
not affects the singlets but flips the up spins. None 
of these A'^ -I- 1 states, are invariant under translation. 
Thus the ground state has a 3(A^ + 1) fold degeneracy, 
generated by the broken symmetries of the Hamiltonian, 
namely the continuous so(3) symmetry and the discrete 
translational symmetry. To find translational invariant 
states, we form the states 1^*^), j = 0, 1, 2 as follows: 
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where lu^ ~ 1, and \(f>i) = T\(f>i^i), where T is the one-site 
translation operator, figure (2). These states are eigen- 
states of T, with eigenvalues l,a; and respectively. 



One can readily show that for i ^ 

from which we find the states \^^) are normahzed in the 
thermodynamic hmit. In this same limit, the correlation 
functions are easy to find. First the magnetization turns 
out to be 
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and the two-point correlation function becomes 
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FIG. 3. The dispersion relation for the holes of model (10). 
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We now ask if it is possible to find some of the exact 
excited states of the Hamiltonian. We will show that 
this is indeed possible, although the states thus found 
are not the low-lying states of the energy spectrum. To 
proceed, we note that the state 
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is an exact energy eigenstate, indeed the maximum en- 
ergy eigenstate with Emax = 6iV J (here we are working 
with h = 2JP2 + Pq and not with (10), in which multi- 
plicative and additive constants have been ignored). To 
see this note that Po\^) =0 and ^2!^^) = 1^^)- The max- 
imality of the energy is found by taking J very large and 
noting the continuity of the energy with respect to J. 
Consider now the one "hole" state of the form 



where the state in the n — th position has been flipped. 
A general one-hole state is of the form 
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It is obvious that Po\n) ~ 0, for any embedding of projec- 
tion Pq on four consecutive sites. To see this, note that 
the total z component of each four consecutive spins is 
either 2 or 1. Also from the explicit form of the spin-2 
states, one can check that P2 does not take a state |n) 
outside the one-hole sector. In fact, a simple calculation 
shows that 
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H\n) = 6J{N - l)|n) + -^^{4: - k){\n + k) + \n - k)). 

fc=i 

Inserting this in the eigenvalue equation H\ip^) = Er\ip^), 
gives 
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FIG. 4. The states in each bulb, are projected to a spin 
3/2 state, to form a new quantum spin chain with spin equal 
to 3/2 and with nearest neighbor interactions. 
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Note that —Emax is the exact ground state energy of 
— iJ, in which case —E^ gives the energy of the low-lying 
excitations or spin waves of — iJ. The dispersion relation 
is plotted in figure (3). Note that if it were not for the 
Pq projector, then a Majumdar- Ghosh state, a singlet, 
would have been the ground state of H . The role of Po 
is to elevate the energy of this state and hence break the 
symmetry. Finally we note that from the above model 
one can obtain an exact ground state for a spin | model 
with nearest neighbor interaction. Let us project three 
spins into a spin | site as shown in figure (4) . The prod- 
uct of two adjacent spin |'s on the lower chain, being ef- 
fectively the product of four spin i's on the upper chain, 
can combine to give only spins 2, 1 and 0. Therefore a 
projector P3 on spin 3 sector, can annihilate the local 
state of any two adjacent spins in the lower chain. This 
means that the state constructed as such, will be an ex- 
act ground state of a Hamiltonian with nearest neighbor 
interaction H = ^^^(-fs)* fc+i- The projector is found in 
the same way as in (7) and (8) and we arrive at 
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-(Sj • Sj+i)^, 



where C = ^N. In this way, we construct a rotationally 
invariant Haniihonian whose exact ground state breaks 
this rotational symmetry. The remarkable point is that 
we can also obtain the explicit form of this state, be- 
yond the construction shown in figure (4) . With this ex- 
plicit expression one can then determine the correlations 
in this state which are no longer restricted to nearest 
neighbor sites. To achieve this, let us first construct a 
Matrix Product (MP) representation for the parent (up- 
per) dimer state in figure (4) . In view of the MP repre- 
sentation [10] of dimerized Majumdar- Ghosh state [11], 
this state has the following representation: 



and 



|0o) = triA. A. B A. A. B 
where 

^+ = l+)(0| + |0)(-|, 

A_ = |-)(o|-|o)(-M, 



I«l,«2,«3, 'i4,«5,i6, ■ • •) 



B_ = 0, 



(18) 



and |-|-), |— ) and |0) are three ortho- normal states. Using 
the Clebsh-Gordon coefficients to decompose the product 
of three states in a bulb of the upper chain in figure (4) 
and noting that B ~ 0, we find the following matrices 
corresponding to the states of a single spin 3/2 site on 
the lower chain: 

4/2 K,. ^ ^{A_B^A^ + A^B^A_), 



^-3/2 = 0, 



leading to 



-^A BA 
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^-3/2 = 0. 



The state of the lower chain is now given by a matrix 
product state, with the above matrices. The ground 
state is a multiplet with total spin N/2, for which we 
have depicted the top state explicitly. This MP repre- 
sentation for spin 3/2 chain has already been obtained 
in [12] by other methods. Having a MP representation, 
it is straightforward to calculate the correlation func- 
tions. The results which agree with that of [12] are that 
(5^) = i, (5^) = {Sy) = 0, and 



((5^'r>=5/4- 



where u := — ^=. 

l + \/3 



{{sn') = {{sn') = 5/4 

With the notation G'"(l,r) 
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G^(l,r) = Aii-iy-' e-^ , G*(l,r) = At{-iy-' 



where 
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(19) 
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(20) 



\n{2 + V3)' ln(l + V3)' 

The correlation functions for other states in the ground 
state multiplet can be determined by the application of 
rotational symmetry. In summary, we have introduced 
a one-parameter family of isotropic spin 1/2 quantum 
chains, whose exact ground states break the translational 
and rotational symmetry of the Hamiltonian. The total- 
ity of degenerate vacua of these models are generated 
by the application of rotation and translation operators. 
Besides the ground state, the maximum energy state and 
some of the lower energy states have also been derived 
in closed analytical form. A related spin 3/2 quantum 
chain with nearest neighbor interactions has also been 
constructed, the exact ground state of which has an MP 
representation. The method presented in this letter can 
be pursued further, to make other spin models whose 
ground state break the rotational symmetry of the Hamil- 
tonian. There is also an interesting un-answered question 
which deserve further investigation: Is there a singlet 
state which is invariant under rotation and has a higher 
energy than the degenerate ground states [13]? V. K. 
would like to thank Laleh Memarzadeh for very valuable 
and constructive discussions. S. A. and S.B. would also 
thank All Rezakhani for valuable discussions. 
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